Hydrodynamical thermotransport relaxation times of semiconductor electrons via 

acoustic phonons 
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We propose a hydrodynamic model to study the thermotransport properties of semiconductor 
electrons. From the semiclassical Boltzmann equation a set of balance equations is built for the 
relevant fields. The electron density, the electron energy density, the electric current density and 
the heat flux density are considered as the basic fields of direct transport and cross effect fluxes. The 
kinetic relaxation times of the production terms are calculated by considering the electron-acoustic 
phonon interaction as the leading microscopic scattering process. To justify the long time thermal- 
ization regime, the Onsager symmetry relations are proved, both on the completely degenerate and 
non degenerate limits. 

PACS numbers: 05.20.Dd, 72.10.Bg, 72.15.Jf, 72.15.Lh 



I. INTRODUCTION 

Hydrodynamical models have been widely used to de- 
scribe thermotransport properties of electrons in semi- 
conductor materials. Complementary contributions of 
charge transport and thermal dissipation fluxes can be 
included in the same model, but due to the complexity 
of the system, and the amount of microscopic processes 
involved, it is mandatory to separate the analysis of the 
microscopic scattering processes of electrons. 

As usual, a generalized hydrodynamic model of semi- 
conductors requires of constructing a hierarchy of balance 
equations for the complete set of basic fields from the 
semiclassical Boltzmann transport equation (BTE)[l)|. 
This hydrodynamic scheme constitutes a closed set of bal- 
ance equations for the relevant variables once the higher 
order moments and the production terms have been com- 
pletely evaluated. 

Depending on the order of the approximation, an 
expansion of the non equilibrium distribution function 
(NEDF) must be used to calculate the higher order mo- 
ments by means of the local variables. For example, in 
degenerate semiconductors, a local Fermi-Dirac distribu- 
tion function is used as a weight function to take in to ac- 
count the non equilibrium regime. Most of the complica- 
tions appear when calculating the production terms cor- 
responding to the non conserved variables. The largest 
difficulty arises as a result of the intrinsic complexity of 
the coUisional kernel of the BTE, enhanced by the mi- 
croscopic semiconductor scattering processes of electrons 
which involve a variety of crystal lattice perturbations. 

Due to all these considerations the closure problem 
is usually solved by including some phcnomcnological 
components, such as the relaxation time approximation, 
or the computational Monte Carlo analysis. The result 
is a hybrid scheme of thermotransport of electrons be- 
tween the kinetic and phenomenological models (or ki- 



netic and numerical Monte Carlo method) 0, . Follow- 
ing these arguments, we can understand why most of 
the hydrodynamical models of thermotransport in semi- 
conductors have some approximations in the region of 
quasi-stationary macroscopic fluxes. In this region, the 
perturbations drive the system not to far from the local 
thermodynamic equilibrium, so it is possible to develop 
a complete kinetic theoretical model considering a spe- 
cific electron-lattice interaction dominating the thermal- 
ization process. 

In this work, a kinetic theory approach is considered to 
solve the semiclassical BTE for the electrons in a semi- 
conductors. The method of moments derives into a set 
of balance equations up to the 8-fields approximation. 
The basic fields of the electron density, electric current 
density, energy density, and heat flux are the local rele- 
vant variables. The closure of the balance equations is 
achieved by using the maximum entropy principle with 
a NEDF, which allows to evaluate the production terms 
without any phenomenological assumption, as is com- 
monly done in the literature P, Q- The electron-acoustic 
phonon interaction is the leading scattering mechanism, 
i.e., the lattice temperature is relatively low. The elec- 
tric applied fleld is absent, however, an internal electric 
fleld is generated because of the inhomogeneties of the 
electron's distribution. The resulting kinetic relaxation 
times are used to justify the linear hydrodynamic model 
by checking the corresponding Onsager symmetry rela- 
tions. This assures that the linear response theory cor- 
responds to the relaxation time approximation. 



II. THE KINETIC MODEL 

We assume that a macroscopic state of the electrons is 
characterized by the eight local relevant variables corre- 
sponding to the number density of the electron's charge. 
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n(x, t), which defines the mass density through p(x, t) = 
mn{x,t), where m is the electron's effective mass. The 
electric current density, Ji{x,t) = —en{x,t) Vi{x,t), 
where Wi(x,i) is the hydrodynamic velocity and e is 
the electrical charge of electrons. The energy density, 
W{x,t), is defined as the average of £ — \'raC^. Here, 
Ck is the peculiar velocity with components Ci (x, k, t) = 
Ci(k) — Vi{:x.,t) and Ci(k) = hki/m is the velocity of elec- 
trons. To complete the set of basic fields, we have de- 
noted by Qi (x, t) the local heat flux density. 

The balance equations for describing the fluxes of elec- 
tron's charge and energy, according to Refs. [l|, are: 



dp 
dt 



m dJi 
e dxi 



= 0, 
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<h_ ^ KbT J, 



The effective local electric field is defined by £i = 
{l/e)d(p/dxi, where cp is the local chemical potential. We 

have difined the constants, ly — 2^/'^K^J'^rTt'^'^T:/h^ , Kb 
the Boltzmann constant, h the Planck constant, and de- 
noted by /„ = Iniifl KbT) the Fermi-Dirac integrals!^. 
Summation over repeated indices has been considered. 
Note that the production terms, rhs of Eqs. (1), contains 
contributions to the cross effect fluxes, which are propor- 
tional to tJq and Tq j respectively. They are necessary 
to compare the constitutive relations of the thermotrans- 
port coefficients under a Onsager symmetry criterion. 

The closure of the set of generalized hydrodynamic evo- 
lution equations is achieved by using a NEDF. This func- 
tion must be an approximate solution of the Boltzmann 
equation coming from the maximum entropy principle, 
which is assumed to be of the form / = f^^{£){l + $). 
The term $ can be written in terms of the basic fields 



where, 



(3) 



ci^KbT^ 



-ij{£)a(k).. 
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(3) 



-/Q(£)a.(k) 
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The functions Ij {£) and Iq {£) are defined through 
the Fermi integrals as follows ^Si], 



and 



5£ 

7 KbT 



25/2 - 21/2/6 



3£ J 



25/| 



21/2/6 



(5) 



(6) 



In this semiclassical hydrodynamical model, the weight 
function f^^{£) is a local Fermi-Dirac distribution func- 
tion. Also, in Ref. |5| we can see that the relaxation times 
have the following general form for the electric current 
density. 



^ = |dkdk'5(k,k')(k'-k) 



(7) 



[fni-fr){i-.f^''')F, 



(3) 



Jk /k' U - Jk' 



and for the heat flux density, 

^ ^H"^ J dkdk'S{k,k')i£'k' ~£k) 
[/k^^(l-/k^^)(l-/k^^)Gf 



(8) 



/^''/k'^''(l-/k'")G, 



FD\ 



Similar expressions have been found for the cross effect 
fluxes. The signs under the summation symbol indicate 
absorption and emmision respectively. 



III. THE RELAXATION TIMES 

To evaluate the kinetic relaxation times, we must spec- 
ify an electron-phonon interaction according to the per- 
turbations in the crystal lattice. In the present case, 
we have considered interactions with acoustic phonons 
as the leading mechanism of the collisional processes of 
electrons. As usual, the rate of electron collisions with 
acoustic phonons is given byQ, 



.(k,k')^| 



DLKbT 



2L^pvl 



S (£ - £') 



(9) 



where Dac is the electron-acoustic phonons deformation 
potential, and Vg represents the velocity of sound. After 
a substitution of Eq. (9) into Eqs. (7-8) and considering 
that the electron energy £ is the integration variable, 
we obtain general expressions for the relaxation times. 
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The corresponding relaxation time for the electric current 
density is, 



^-1 ^ 7 {27:fDl,m^ 



^ 2 h^pv1{KBT)^ 

jlj{£) [f°{£)f{l-f°{8))8'd8 



(10) 



A similar procedure gives rise to all the other relax- 
ation time expressions. The following are the expressions 
for l/rg, 1/tjq, and 1/tqj. For the heat flux we have, 



, 5 (2^)^Z?2^mi 

= - 



2 hipvl{KBT)i 
jlQ{E) [f^^{S)\'{l~f^^{E))E'dE 



(11) 



and the crossed effect relaxation times, tjq and tqj, are 
given by, 



.1 _ 5 {2TTfDl^mi 



(12) 



Iq{E) [f^^{£)Y{l-f^^{£)) E'dE 



7 {2nfDlmi 
2 h^pv1{KBT)i 

Ij{E) [f^'^iE)]' 



(13) 



I- f''°{E))E^dE. 



By incorporating this expressions of the relaxation times 
in the hydrodynamical set of transport equations we ob- 
tain a completely closed model containing no free or phe- 
nomenological parameters. It is very important to note 
that though the structure of the transport equations is 
given, and does not depend on a microscopic model for 
the interaction, the specific values of the relaxation times 
depend closely on the model of the electron- lattice inter- 
action. 



IV. ONSAGER SYMMETRY 

As is well known, the Onsager symmetry principle 
projects microscopic reversibility of the interactions onto 
a macroscopic scale through the transport symmetry re- 
lations. In this section we will show that, in the linear 
response regime, the semiconductor transport of charge, 
under internal electric fields variations and thermal gra- 
dients, can be characterized by a set of constitutive rela- 
tions with thermoelectric transport coefficients display- 
ing symmetry on the cross effect fluxes. This point de- 
serves some discussion because the microscopic interac- 
tions of the semiconductor electrons are regulated by a 
process of creation and annihilation of phonons. So there 



is no guarantee of microscopic reversibity. However, be- 
cause of the elastic collision between electrons and acous- 
tic phonons, our result reflects the existence of a kind of 
a more general detailed balance principle of collisions. 

We show the validity of the Onsager symmetry rela- 
tions, at the end of the thermalization process and both 
in the completely degenerate as well non degenerate lim- 
its, by following a general expression that must be satis- 
fied (see Ref. Ill): 



Tj ( 5/4-3/2;^ ) = -TQ 



5h-7Ie—]. (14) 



First, we are considering the kinetic relaxation times in 
the complete degenerate limit. In this case, integration 
with respect to energy, and assuming the Fermi-Dirac 
distribution function satisfying f^^{E) (l — f^^{E)) = 
d (E — If), derives into the following expressions for the 
relaxation times: 



5y T, 

7KbT^'*: 



7 h_ 
2 25/| - 21/2/6 ' 
5 -^4 - 5^/2 



(15) 



2 ^ 25/|- 21/2/6 /^sT' 
2 25/| - 21/2/6 ' 



h - 



25/?- 21/2/6 /^sT' 



Here, we have introduced the parameter v, which de- 
pends on the chemical potential of electrons and the local 
temperature. 



/ 2nmD. 



\ KbT 



KbT KbT 



(16) 



A direct substitution of these expressions into the On- 
sager conditions show us that they are well satisfied. 
On the non degenerate case, the Onsager symmetry re- 
lation is well satisfied too. If we consider this limit, 
f^^{E) <C 1 and it becomes a Maxwellian distribution 
function. Under these conditions the relaxation times 
can be expressed as follows. 
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here we have defined the parameters, 

{2T:Dacfmi ( 2Lp 
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Once again, the Onsager condition is well satisfied in 
this case. Next, what we can see is if the full integral 
expressions for the kinetic relaxation times satisfy the 
symmetry condition. However, this requires a computa- 
tional analysis and a correct interpretation of the numer- 
ical uncertainties. We can speculate that the Onsager 
symmetry is going to be well satisfied in all the interval 
of the temperature varations, because of the regime in 
which we are working (the last stage of thermalization) 
and the assumption of the microscopic detailed balance 
under elastic electron-phonon collision being satisfied. 

The general condition of symmetry, in Eq. (14), also 
appears in Refs. [Ell], among others, with respect to the 
Onsager symmetry relations of the thermotransport coef- 
ficients. In particular in Ref. it is shown that for ther- 
motransport in polar semiconductors, where the main 
microscopic interaction mechanism is the electron-polar 
optical phonons, the relation specifies that only some val- 
ues of Lp are satisfying the symmetry condition, i.e., only 
certain values are allowed for the chemical potential in 
comparison to the optical phonon energy. Though this 
condition seems a very restrictive one, we must note that 
the model corresponds only to the behavior of electron 
population. The phonons become the reservoir of lat- 
tice perturbations, so they are not taken into account in 
order to preserve microscopic reversibility:?]. Hence the 
Onsager symmetry relation depends strongly on the ther- 
malization regime where a different microscopic electron- 
lattice interactions may be dominating the relaxation 
processes. 

V. DISCUSSION AND CONCLUSIONS 

We have developed a kinetic theoretical approach to 
describe a generalized hydrodynamic behavior of semi- 
conductor electrons at the end of the thermalization pro- 
cess. The semiclassical Boltzmann transport equation 
have been solved up to the eight moments method ap- 
proximation in such a way that the heat flux density is 
considered as a relevant variable in the resulting hydro- 
dynamical model. The closure of the balance equations 
is achieved by evaluating the higher order moments, as 
well as the production terms, through a non equilibrium 
distribution function according to the maximum entropy 
principle as has been described in literature. In case that 
the electron-acoustic phonon scattering is the leading mi- 



croscopic interaction, the production terms can be ex- 
pressed by means of a set of kinetic relaxation times. 

As a consequence of the theoretical model developed 
in this paper, we have shown analytically that condition 
(14) is satisfied in both degenerate and non degenerate 
limits. This means that Onsager symmetry relations are 
fulfilled. Even more, and given that for the electron- 
acoustic phonon scattering the microscopic reversibilty 
is valid, we believe that the general expressions for the 
relaxation times, Eqs. (10)-(13), will give that Onsager 
symmetry relations are also satisfied under more general 
conditions. 

In general, the electron transport properties in semi- 
conductors depend on the scattering processes that elec- 
trons suffer. In addition, these transport processes can be 
expressed in terms of the relaxation times for each scat- 
tering mechanism. Relaxation times are usually obtained 
by Monte Carlo Simulations^]. For the case of electron- 
acoustic phonon interaction we have found closed gen- 
eral expressions, Eqs. (10)-(13), which_ together with 
the general expressions found in Ref. for electron- 
optical phonon interaction provide us with a theoreti- 
cal framework to describe the thermotransport processes 
of electrons in semiconductors under more general con- 
ditions with no free or phenomenological parameters. 
A complete scheme will be obtained when we include 
electron-intervalley phonons and electron-impurity scat- 
tering processes. The case of piezoelectric phonon per- 
turbations is rather direct and follows the same argument 
as acoustic phonon perturbations. A detailed analysis of 
the relaxation times as function of lattice temperature 
and electron concentration will give us useful informa- 
tion of the main relaxation channel under the imposed 
conditions 0, [13] ■ For example, in Ref. [§] it is shown 
that, for a silicon device, the energy relaxation time 
(which is the characteristic time in which the electron 
energy relaxes to the lattice) is substantially larger than 
the relaxation times for momentum and heat flux. The 
situation could be very different for Bulk GaAs. For high 
temperatures (300 K) and low electron concentration (of 
the order or less than 10^^ cni^'') the main collision 
mechanism is electron-acoustic interactions. Whereas, 
for low temperatures and high electron concentration 
where screening effect should be taken into account, the 
electron-optical phonon interaction should be the main 
ralaxation mechanism Jj. 
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